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ABSTRACT: Tight focusing of double ring shaped radially polarized multi gaussian beam with cosine phase plate is 
studied numerically based on vector diffraction theory. For higher-order, radially polarized mode beams as well as a 
fundamental mode (R-TEMo1*) beam, the strong longitudinal component forms a sharper spot at the focal point 
under a high-NA focusing condition. In particular, double-ring-shaped radially polarized mode (R-TEM11*) beams 
can effectively reduce the focal spot size because of destructive interference between the inner and the outer rings. 
The mathematical expressions for the focused fields are derived. Simulation results show that the focused fields and 
phase distributions at focus are largely influenced by both the cosine parameter and beam diameter of the incident 
beams. Moreover, focal spot with flat-topped, two focal spots patterns can be flexibly achieved by carefully choosing 
the cosine parameter (C) and the ratios of the pupil diameter to the beam diameter (8) which confirms the potential 


of such beams in wide applications, such as optical tweezers, lithography, and material processing. 
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1. INTRODUCTION 

Tight focusing of radially polarized light is a very 
topical area of research [1-3]. It has been shown that, 
in many circumstances, focusing of radially polarized 
light results in a smaller focal spot than focusing 
linearly polarized light. Applications are in microscopy, 
lithography, metrology, nonlinear optics, optical 
manipulation, and optical micromachining. Gaussian 
beams are very similar to researchers and have found 
wide applications in many optical systems; however, 
this kind of beam is only one solution of the Helmholtz 
equation. In fact, Gaussian beam is one special case of 
Hermite-sinusoidal-Gaussian beams that are more 
general solution and were introduced by Casperson 
et.al. [4-6]. Hyperolic-cosine-Gaussian (cosh-Gaussian) 
beams are regarded as the special case of Hermite- 
sinusoidal Gaussian beams, and Gaussian beams are 
still base cases of cosh-Gaussian beams. The cosh- 
Gaussian beams have attracted much attention because 
their profiles can be altered by choosing the suitable 
decentered parameters in cosh parts [7-9], and 
propagation and focusing properties of cosh-Gaussian 
beams have become object of many works [10]. Zhang 
et al proposed the virtual sources for a cosh-Gaussian 
beam [11]. Zhang et al proposed the virtual sources for 
a multi Gaussian beam [12]. However, most of the 
studies on the radially polarized beam and _ its 
applications have concerned the so-called doughnut- 
shaped beam, which is represented as (TMoi mode) 
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[13, 14]. Recently, Kozawa et al report that the radially 
polarized higher-order Laguerre-Gaussian beam (LG) 
has higher a fundamental radially polarized mode 
performance than the low-order one and can generate 
the focal spot beyond the Abbe's diffraction limit 0.5A. 
This kind of higher-order LG beam has concentric multi 
ring-shaped intensity pattern with m phase shift 
between the rings[15-17]. 

In the optical trapping system, it is usually deemed 
that the forces exerted on the particles in light field 
include two kinds of forces; one is the gradient force, 
which is proportional to the intensity gradient, and the 
other is the scattering force, which is proportional to 
the optical intensity [18, 19]. Therefore, the tunable 
focal shift predicts that the position of optical trap may 
be controllable. In fact, tracing the movement of the 
point of absolute maximum intensity along the axis has 
become the subject of a number of publications for 
several decades [20-23]. It was found that the point of 
absolute maximum intensity does not coincide with the 
geometrical focus but shifts along the axis towards the 
aperture of the optical system [22-24]. More 
interestingly, the focal shift may be in continuous in 
certain optical focusing systems. It was found that the 
focal shift may be accompanied by an effective 
permutation of the focal point, and this effect is 
referred to as focal switch [25,26]. Focal shifts have 
taken an important part in design and optimization of 
various optical systems and extremely important for 
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image formation in defocused planes, for increasing 
the depth of focus and in automatic focusing [27]. More 
interestingly, the focal shift may be in continuous in 
certain optical focusing systems. It was found that the 
focal shift may be accompanied by an effective 
permutation of the focal point, and this effect is 


the focusing properties of the double ring shaped 
radially polarized multi gaussian beam with radial 
cosine phase plate are not studied. In this paper the 
focal shift of the double ring shaped radially polarized 
multi gaussian beam with radial cosine phase plate is 
investigated theoretically by vector diffraction theory. 


referred to as focal switch [28,29]. In our knowledge, 


2. VECTOR DIFFRACTION THEORY OF DOUBLE RING SHAPED RADIALLY POLARIZED MULTI 
GAUSSIAN BEAM 

A schematic configuration of the proposed method is shown in Fig.1. An incident double ring shaped radially 
polarized multi gaussian beam transmits through a radial cosine phase plate and is subsequently focused by a 
high NA lens. The analysis was performed on the basis of Richards and Wolf's vectorial diffraction method [30] 
widely used for high-NA lens system at arbitrary incident polarization. In the case of the incident polarization, 
adopting the cylindrical coordinates r,z,@ and the notations of Ref. [31], the electric field E(r,z,) in the vicinity 
of the focal region can be written as 


E(*z)= 


are the amplitudes of the two orthogonal components and @,and @. 


E,é,+Ezéz > 


Where E, and E. are their 


corresponding unit vectors. The two orthogonal components of the electric field is given as 
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Where @ = arcsin( NA) / 7 the maximal angle is determined by the numerical aperture of the objective lens, 


and n is the index of refraction between the lens and the sample. k = 277/ J is the wave number and Jn(x) is 
the Bessel function of the first kind with order n. r and z are the radial and z coordinates of observation point in 
focal region, respectively. @ describes the pupil apodization function. The intensity in the focal region is 
proportional to the square modulus of Eq. (1). In the system investigated in this chapter, the wavefront phase 
distribution is radial cosine function distribution, and can be written as [32] 


tan(0) | 
tan(c) | 
where C is the frequency parameter in cosine part of the wavefront phase distribution, C denotes the radial 
change frequency of the phase and parameter. The reason for choosing this kind of radial cosine phase 
wavefront is that it is very simple and easy to carry out, for example this kind of phase distribution can be 
implemented by phase spatial light modulator or by pure phase plate manufactured by the lithographic method 
conveniently. The electric field in the focal region can be written in the same form as Eq. (1); however, the three 
orthogonal components Er, and Ez are different and should be expressed as 
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The total intensity distribution in the focal region is proportional to the square modulus of Eq. (1), the radial 
polarized component and longitudinal polarized component; azimuthally polarized component can be 
calculated according to Eqs. (5) and (6) respectively. The 7(@) describes the double ring shaped radially 


polarized multi gaussian beam, this function is given by [12] 
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Here, 0 is the converging semi-angle. We denote the maximum converging semi-angle as Omax which is related 
to objective numerical aperture by O@max= arcsin(NA). 80 is an angle which, along with integer m, determines the 
shape of the modulation function.@ois usually chosen to be slightly smaller than O@max. 9c determines the radial 
position translation of P(@). Here we take 0c =@max/2. wo is the waist width of single Gaussian beam which is 
calculated by the following formula 


Omax 


weft ; eon)! 


Equation (7) describes an object beam. It is suitable to convert Eq. (8) from a function of angle into a function 
of radial polar coordinate. One may substitute 0 with arcsin(r/f), where f is the focal distance of the objective. 
In Eq. (7), the factor (0/0) measures the most of light energy is located on the annular edge of the pupil. 
Increasing the integer m concentrates more energy into the annular edge area in which the converging semi- 
angle is more than 60. The sum of (2N + 1) spatially equally spaced Gaussian beams ensures that amplitude of 
the constructed annular multi-Gaussian beam decreases suddenly, when reaching the outer edge of the pupil. 
T(0@) denotes the apodization function of the LGi:1 beam, which is given as [15-17], 


2 2 
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Here Lp is the generalized Laguerre polynomial and 8 is the ratio of the pupil radius to the incident beam 
radius in front of the focusing lens. Note that B should be greater than 1 because the outer ring of the R-TEM 11 
beam will be completely truncated by the pupil if Bs1. In this study the NA is assumed to be 1.2, with n=1.33.A 


is determined such that each truncated beam has the same optical power. 


3. RESULTS AND DISCUSSION 

Without loss of validity and generality, the 
focusing properties of the double ring shaped 
radially polarized multi Gaussian beam with radial 
cosine phase wave front are calculated. We perform 
the integration of Eq. (1) numerically using 
parameters A = 1, wo=0.125, NA = 0.9 and m=1. Here, 
for simplicity, we assume that the refractive index n 
=1 and A=1. For all calculation in the length unit is 
normalized to A and the energy density is normalized 
to unity. In order to understand the focusing 
properties of the double ring shaped radially 
polarized multi Gaussian beam extensively, firstly, 
the focusing of double ring shaped radially polarized 
multi Gaussian beam with high NA lens system is 
investigated without wave front phase modulation. 
The intensity distribution of a radially polarized 
multi Gaussian incident beam for the high NA lens is 
shown in Fig.2 (a-c). It is observed from the Fig.2a 
shows the intensity profile of the radial, longitudinal 
and total electric field components of the optical field 
at focus. It is evident that the intensity of the 
longitudinal component (blue line) is higher than the 
radial component (red line). From the Fig.2a, we 
measured the FWHM of the generated focal spot size 
is 1.24 and focal depth 2A which is shown in Fig.2 b. 
Such a focal spot segment is useful for high refractive 
index particle trapping. The Fig.3 (a, b & c) shows the 
focal segment generated for 6 =1, for cosine 
parameter values of C=0.25, C=0.5, and C=1. It can be 
seen from the Fig.3 (a) shows that the generated 
focal segment is a shifted focal spot. From the Fig.3 
(a) the maximum intensity position shifted to -0.9A, 
However when we increased cosine parameter (C) 


values for C=0.5 and C=1 the maximum intensity 
position shifted to -1A , and -2A_ respectively is 
shown in fig.3 (b & c). We observed that the 
generated focal spot in focal region shifts towards 
the optical aperture along optical axis on increasing 
cosine parameter C, namely focal shift phenomenon 
occurs.Figure.4 shows the intensity distributions of 
double ring shaped radially polarized multi gaussian 
beam without cosine plate for B=1.26. As shown in 
Fig.4(a), we measured FWHM of the splited focal 
spot is 0.8 A and corresponding radial component 
percentage is 70% of the total intensity and the 
FWHM of the sub wavelength focal spot size is about 
0.251, and focal depth of 2.8 A as shown in Fig. 4b, 
axially separated distance is 2.8 A is shown in fig.4 
(c). Such type of focal spot used for multiple optical 
trapping applications. The Fig.5 (a, b & c) shows the 
focal segment generated for B =1.26, for cosine 
parameter values of C=0.25, C=0.5, and C=1. It can be 
seen from Fig.5 (a) shows that the generated focal 
segment is a two focal spots. From the Fig.5(a-c) we 
observed that the generated focal spots in focal 
region shifts towards the optical aperture along 
optical axis on increasing cosine parameter C. 
However when we change the £ value from 1 to 1.26 
the single focal spot is changed to two focal spots 
segment. The Fig.6 shows the focal structure 
generated by double ring shaped multi gaussian 
beam tightly focused with high NA lens. It is 
observed from Fig 6(a), when B= 1.08 and the cosine 
parameter value C=0, the generated focal generated 
focal structure is a flattop profile having FWHM of 
1.6A and focal depth of 3.8A.The 2D _ intensity 
calculated at the focus illustrates that the focal 
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structure is dominated by the longitudinal 
component and the radial intensity remain as a 
intensity as 80% of the total intensity. Now, the 
beam parameter f is chosen as 1.08 to investigate its 
effect on focal shift. The intensity distributions for 
(G=1.08) and different C are illustrated in Fig.7. It can 
be seen from this figure that the beam parameter C 
affects focal intensity distribution very considerably. 
On increasing phase parameter C, one focal peak 
evolves into one dark focus that is defined as those 
whose intensity is lower than that around it, namely, 


cosine phase plate 


Double ring shaped 
radially polarized Multi 
Gaussian beam 


6... aX 
7 \ _P(t.9,z) 
a 
. Focal Plane 
<————_ @ 


local focus intensity minimum. There is one optical 
intensity maximums along optical axis. On increasing 
phase parameter C, whole focal pattern shifts along 
optical axis near from the optical aperture, one 
intensity peak shrinks so that there is only one 
optical intensity maximum. The dependence of focal 
shift on C is also given in Fig.7(a-c). If this focal peak 
is used to construct one optical trap, cosine 
parameter C may be employed to adjust trap 
position, in turn, can transport micro particles. 


Ny 


— 


High NA Lens 


Fig.1. Focusing of a double ring shaped radially polarized multi Gaussian beam with radial cosine phase 
plate by high NA lens 
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Fig.2. (a) 2D intensity distribution at z =0A. (b) On axial intensity at r= 0A (c) 3D intensity distributions 
in r-z plane incident R-TEM 11* beam for B=1 without wave front phase modulation. 
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Normalized total intensity distributions in r-z plane incident R-TEM 11* beam under the condition 


of B=1, (a) C= 0.25, (b) C = 0.5(c ) C = 1 respectively. 
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Fig.4. (a) 2D intensity distribution at z =0A. (b) On axial intensity at r= 0A (c) 3D intensity distributions 
in r-z plane incident R-TEM 11* beam under the condition of B=1.26 without wave front phase 
modulation. 
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Fig.5. Normalized total intensity distributions in r-z plane incident R-TEM 11* beam under the condition 
of B=1.26, (a) C = 0.25, (b) C= 0.5(c ) C= 1 respectively. 
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Fig.6. (a) 2D intensity distribution at z =0A. (b) On axial intensity at r= 0A (c) 3D intensity distributions 
in r-z plane incident R-TEM 11* beam under the condition of B=1.08 without wave front phase 
modulation. 
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Fig.7. Normalized total intensity distributions in r-z plane incident R-TEM 11* beam under the condition 
of B=1.08, (a) C= 0.25, (b) C=0.5(c ) C = 1 respectively. 


4. CONCLUSION 

In conclusion, the focusing properties of the 
double ring shaped radially polarized multi Gaussian 
beam with radial cosine wavefront phase 
distribution are studied based on vector diffraction 
theory. The intensity distribution in the focal region 
is investigated in detail by numerical calculations for 
different cosine parameter values (C) and truncation 
parameter (8). The simulation results show that the 
intensity distribution in focal region of the double 
ring shaped radially polarized multi gaussian beam 
can be adjusted considerably by the cosine 
parameter C. Focus can shift along optical axis on 
increasing C, and focal pattern changes are observed. 
Moreover, by altering frequency parameter of the 
phase plate will change the energy distributions of 
maximum intensity peak and other small intensity 
peaks appeared and novel focal patterns also evolve 
considerably. The focusing properties technique may 
find wide applications in some domains such as 
optical tweezers, optical storage, particle 
manipulation, multifocal microscopy, laser printing. 
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